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Polymer quantization and the saddle point approximation of partition functions 

Hugo A. Morales- Tecotl, 1 ’^ Daniel H. Orozco-Borunda , 1 and Saeed Rastgoo 1 [] 

1 Departamento de Fisica, Universidad Autonoma Metropolitana - Iztapalapa 
San Rafael Atlixco 186, Mexico D.F. 09340, Mexico 
(Dated: November 11, 2015) 

The saddle point approximation of the path integral partition functions is an important way 
of deriving the thermodynamical properties of black holes. However, there are certain black hole 
models and some mathematically analog mechanical models for which this method cannot be applied 
directly. This is due to the fact that their action evaluated on a classical solution is not finite and 
its first variation does not vanish for all consistent boundary conditions. These problems can be 
dealt with by adding a counterterm to the classical action, which is a solution of the corresponding 
Hamilton-Jacobi equation. 

In this work we study the effects of polymer quantization on a mechanical model presenting 
the aforementioned difficulties and contrast it with the above counterterm method. This type of 
quantization for mechanical models is motivated by the loop quantization of gravity which is known 
to play a role in the thermodynamics of black hole systems. 

The model we consider is a nonrelativistic particle in an inverse square potential, and analyze 
two polarizations of the polymer quantization in which either the position or the momentum is 
discrete. In the former case, Thiemann’s regularization is applied to represent the inverse power 
potential but we still need to incorporate the Hamilton-Jacobi counterterm which is now modified 
by polymer corrections. In the latter, momentum discrete case however, such regularization could 
not be implemented. Yet, remarkably, owing to the fact that the position is bounded, we do not 
need a Hamilton-Jacobi counterterm in order to have a well-defined saddle point approximation. 

Further developments and extensions are commented upon in the discussion. 
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I. INTRODUCTION 


Two major open problems in theoretical physics regard the nature of spacetime: on one hand there is the issue of 
singularities, beyond which classical general relativity cannot be continued, and, on the other, one finds the divergent 
high energy behavior of field theories. A quantum theory of gravity is expected to have a bearing on both of these 
problems. For instance, loop quantum gravity jlm l has been shown indeed to replace the big bang of general relativistic 
homogeneous cosmological models by a bounce^ 7| and quantum field theory in such a scenario jH is rather different 
from usual fixed background field theory. Moreover, it has also provided a specific account for black hole entropy Q. 

Now to investigate the behavior of some simple systems under this approach it is possible to use polymer quantum 
mechanics, a finite number of degrees of freedom scheme including some of the loop quantum gravity techniques 3- 
This simplified approach has been applied to some systems to contrast their behavior with their usual Schrodinger 
quantization and its relation with the latter either as a continuum limit 33 or as a low energy approximation 
|13| . Furthermore, for the case of fields on a fixed background this technique has been applied to each of its infinite 
modes as a candidate to explore their high e nerg y behavior fbil . [l5 ] . Even higher order derivative models have been 
given consideration along these lines recently [l6j . With the exception of [17] that advanced a path integral version of 
polymer quantum mechanics, most work on these lines adopted a Hamiltonian scheme. This changed recently: a path 
integral approach was considered [18|-I20] in order to provide a more detailed link between loop quantum cosmology 
and the covariant spin foam models [5J, and a polymer path integral in quantum field theory and its relation to Lorentz 
invariance has also been studied in |2l| . Also the Feynman formula for other mechanical examples has been worked 
out [22], and explicit polymer propagators have been obtained (HI]. An interesting aspect of this so-called Feynman 
approach is that semiclassical approximations can be at hand to investigate important gravitating systems through 
the saddle point approximation of its path integral description. In fact, the Euclidean path integral is specially useful 
in studying the thermodynamics of systems such as black holes since there it is interpreted as the partition function 
of the system in a canonical ensemble. 

Let us consider the semiclassical approximation for a system with the Euclidean path integral 


Z = 




( 1 . 1 ) 


in which <pj are the fields of the theory and S E is its Euclidean action. Given that one can expand the action around 
the classical solutions cf>j | d as 

Se [<Ail c / + Hi] = Se |>|J + SS E [^\ cl + Hj] + \s 2 S E [^\ cl +«%] + • • • , (1.2) 

one can then substitute this into ED; by keeping only up to the quadratic term, one gets 

z ~ ex p (~1 Se nj) / (-^^ 2se [4>i\d+Hi]j, (!-3) 

which is the saddle point approximation to the model that gives us access to the semiclassical information about the 
system. 

More precisely, a saddle point approximation (II.HI) of ED is possible only if the following conditions are met: 

1. The variational principle is well-defined: S E must be functionally differentiable for all the variations of the fields 
and compatible with the boundary and falloff conditions of the fields, so that any boundary term must vanish 
by virtue of these conditions; thus we can write 5S E = f d n x^^-5(j)j. This is needed so that one is able to 

(a) expand S E as in (11.21) around an arbitrary configuration cf )q using the functional derivatives, 

(b) find the classical solutions by setting SS E = 0, 

(c) perform the expansion (11.21) specifically around the classical solutions 4>j\ c i- 

2. Given condition 1, then S E on classical solutions must remain finite, i.e., *S , -e| iSSe _ 0 < oo. 

3. Given conditions 1 and 2, the Gaussian integral jYlj ex P %& 2 S E [ 4>j\ c i + Hj\) also must remain finite. 

4. Since there is a minus sign in the exponent of ED, we should also have S 2 S E \<f>j\ l + S4>j\ > 0 so that the 
classical solutions give the dominant contribution to the saddle point approximation. 
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The physical reason behind the above conditions is that for a semiclassical approximation, the most important 
contribution comes from the classical solutions, and thus everything is expanded around such a trajectory. The rest of 
the terms become less and less important, and thus we only need to keep the perturbative terms up to the quadratic 
term, which gives us the nonclassical contributions. 

There are, however, some important systems for which these conditions are not met and thus access to the semi¬ 
classical approximation via the saddle point method is not possible [24! [29l |. Among these systems are a class of two- 
dimensional (2D) dilatonic gravitational systems, including models like CGHS (Callan-Giddings-Harvey-Strominger 
model) [3Q| and 3+1 spherically symmetric, which have black hole solutions and many other related interesting 
properties. This class can be described by the generic action [3lj | 

5 = -—V/ <Pxy/=jj[$R-U($)V a $V a $- f dxy/qQK. (1.4) 

lb7rG 2 J M ottG 2 Jqm 

Here $ is the dilaton held, and £/(<!>) and H(<h) are model-dependent functions of the dilaton held. The latter is called 
the dilaton potential. The boundary term is the equivalent of the Gibbons-Hawking-York term [32] in this theory in 
which K is the trace of the extrinsic curvature of the boundary manifold. 

To see the problem with saddle point approximation, we consider the on-shell variation of the action j3l| 

SS\, = [ dxy/q ( E ab Sq ab + . (1.5) 

JdM 

It can be shown that this variation does not vanish in some of the models of this generic class. This is basically due 
to the fact that in these models, the coefficients of the variations of the spatial metric q a b and the dilaton $ in the 
above expression, diverge more rapidly than the variations themselves fall off. This might look odd since the presence 
of the Gibbons-Hawking-York (GHY) term is supposed to guarantee the well-posedness of the variational principle in 
these kinds of theories. Notice, however, that the presence of the GHY terms is to let us avoid prescribing Neumann 
boundary conditions for the metric, i.e. it cancels all the variations S(d a gbc) that come from the bulk term. It does 
not guarantee that the Dirichlet boundary conditions on the dilaton field lead to a well-posed variational principle. 
Also it may not be helpful in dealing with issues emanating from falloff conditions of the dilaton field. It is these 
types of boundary and falloff conditions that contribute to the problem here. 

In some of the submodels of this class, even if the variational principle is well defined, the on-shell value of the 
action diverges, especially due to the falloff conditions on the dilaton field and its value on the boundary. So the 
saddle point approximation collapses for these types of models. 

Since as mentioned above, one important use of saddle point approximation is to study the thermodynamics of 
black holes, not being able to make such an approximation for this class of 2D models is a significant shortcoming. 
Luckily there is a common and rather generic method of fixing this problem that consists of the addition of a boundary 
counterterm to Se, which not only makes the action functionally differentiable for all boundary and falloff conditions, 
but also renders the on-shell value of the action finite. It turns out that this boundary counterterm is a solution to the 
Hamilton-Jacobi equation of the system [31]. The action that is the sum of (II. Jl) and the Hamilton-Jacobi boundary 
term is called the improved action. An interesting observation is that the improved action actually gives the correct 
thermodynamics for the system 0. 

In this paper, we investigate the behavior of the polymer path integral quantization of a mechanical model that in 
the usual path integral quantization suffers from the aforementioned ill-defined semiclassical approximation. Whether 
this fixes the problem, and whether we need a Hamilton-Jacobi counterterm, it is relevant to know how the polymer 
approach changes the partition function for such systems and, perhaps, the thermodynamics of some of the 2D black 
holes. The present work is a first step in this direction, and we study a simpler analog model of the aforementioned 
2D class, namely, a particle in an inverse square potential, that has the same technical problems within the usual 
nonpolymer path integral quantization. Let us notice that the Hamiltonian polymer quantization of a particle subject 
to a Coulomb |.33l | and an inverse squared j34j potential, respectively, have been considered previously. Hence, our 
work complements these studies. 

The paper is organized as follows. In Sec. El we introduce the analog model and describe the saddle point issue 
for this case. Then we recall its solution through the application of the Hamilton-Jacobi method. In Sec. Mil we 
polymerize this analog model in a polarization where the position q is discrete. We will see that in this polarization, one 
still needs to add a counterterm to the action to get a new well-defined action suitable for saddle point approximation. 
However, we show that this counterterm and the bulk action are both modified by polymer quantization. In Sec. 
m we propose an effective potential for the semiclassical version of the analog model in a polarization with discrete 
momentum p. We then describe how this effective potential does not require the addition of any counterterm to the 
action. Finally we summarize and make our concluding remarks in Sec. [Vj Details of the calculations are given in the 
appendixes at the end of the paper. 
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II. ANALOG MECHANICAL MODEL AND ITS IMPROVED ACTION 


There are several types of simpler analog models that exhibit the aforementioned ill-defined semiclassical approxi¬ 
mation. One such class of models corresponds to single particle systems in half-binding potentials W j35|. We choose 
a simple system in this class with an inverse square potential with the Newtonian action (the subscripts N and E 
stand for Newtonian and Euclidean, respectively) 

sn= C7 dtLN= Co°° dt ^{^) 2 - w(q) )' <2i> 

and corresponding Euclidean action (by a Wick rotation t —> r = it), Sn = iS E , 


Se 



dr Le 


rTf —i' 

/ 

J Ti—iO 


dr 





The Newtonian equation of motion is 


d 2 8W ( 2 k\ 

Tn & q = -dj-\Tf)' 

which under a Wick rotation becomes 

d 2 8W f 2k\ 

m d^ q = -^{=-7)- 


( 2 . 2 ) 


(2.3) 


(2.4) 


Let us see the problems of the saddle point approximation for this analog model. First, we consider the variation of 
the Euclidean action 


SSe 


OO 




d 8 Le 
drd£q 



Sq. 


(2.5) 


If the boundary term does not vanish for all the variations of Sq compatible with the boundary and falloff conditions, 
then the action is not functionally differentiable. It turns out that this is the case. At r = iti = iO, since the value of 
q is finite, we have Sq = 0. However, at r = itf = ioo, since q —> oo, the condition Sq | ^ 0 is also allowed; i.e., any 

two trajectories do not necessarily coincide at infinity, and yet they both tend to infinite values. Thus the action is 
not functionally differentiable with this boundary condition. 

The common way to overcome this issue is to add a boundary term — B to the action that cancels out the boundary 
term present in (12.51) . Clearly the variation of the boundary term should obey SB = Sq = pSq. It just happens 

dr ^ 

that Hamilton’s principal function G has exactly this property. This is because we have 


SG(q,t) 


dG X A 
— Sq = pdq 

oq 



Sq. 


( 2 . 6 ) 


Hence we conclude that by adding 


~B = -G 

to the action, it becomes functionally differentiable. Such an action is called the “improved” action, 


T[q\ =S E - G(q,t)\l f = dt + ^ - G(q,t)\ t 0 t . 


Clearly, the variation of 
SS E = 


yields 


/ d 8 L e 8L e ) j. 

*l-*5jc; + -srr- 


8L e 8G\ 

WG~G 5q 


J 0 


j ( d 8 L e 8L e j x 

dT h rWq + ^) 5q - 


(2.7) 


( 2 . 8 ) 


(2.9) 


Next, we consider the Euclidean action itself. We would like to show that even if the action is functionally differentiable, 
the value of the Euclidean action on classical solutions is not finite, and that the addition of G to the action, makes 
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it finite. Assuming for the moment that Se is functionally differentiable, using the Leibniz rule on the kinetic term 
and then computing the action on the classical solutions, we can write 


Se\ = 

I 9c! 


m dq 


r=iO 


pic 

/ 

JiO 


. m d 2 q 

dT I _ ' 2 ,d * 3 


k 

-)- 2 

cl <lcl 


( 2 . 10 ) 


Using (12.41) and the form of the classical solution q c i[r), the integral turns out to be finite and of order 7r. On the 
other hand, it can clearly be seen that the boundary term above (or its Newtonian counterpart) diverges, since at 
r — y oo (or t — > oo) the system behaves as a free particle, ^y| c ; (or ^|| c i) is constant, and q c i —> oo. Thus the action 
evaluated on classical solutions is not finite, 


I 


—> oo. 


( 2 . 11 ) 


These results show that the direct saddle point approximation cannot be performed for this system. 

Now let us see what is the form of G and how it makes the action, evaluated on classical solutions, finite. The 
variable t is used from now on, both as Euclidean or Newtonian time, depending on the context. Since G in (12.81) is 
a solution of the Hamilton-Jacobi equation, we have 


H[q, 


dG 

dq 


dG 

~dt 


= 0. 


To find the explicit form of G we need to solve (12.121) . Considering the Euclidean Hamiltonian 

- 2 k 


and noting that p = ^ , the Eq. (12.121) becomes 


2m \ dq 


H = 


dG 


_ 

2 to q 2 ’ 


dG 

~dt 


= 0. 


( 2 . 12 ) 


(2.13) 


(2.14) 


Since the Hamiltonian does not depend explicitly on t, Eq. (12.121) implies that the derivative dG /dt must be a constant 
in time and is actually the negative of the energy of the system H = E. So one can make an additive separation of 
the variables in G as 


which turns (|2.141) into 


G(q, E, t) = fl(q, E) — Et, 

(2.15) 

1 

fdn \ 2 k 

(2.16) 

2m ^ 


W 

\j2m\J Eq 2 + k 

d 

(2.17) 


This can be integrated to yield f2, which together with (12.151) leads to 


G(g, t) = Cq — Et + V2y/rny/ Eq 2 + k — \f2\fmk tan 


-l 


Eq 2 


(2.18) 


Since at q — > oo the potential is zero, the energy takes the form E = ^ m ^ of a free particle (which is equal to p 2 /2m 
asymptotically), thus Hamilton’s principal function behaves asymptotically like 


mq 

2 t ‘ 


(2.19) 


Noticing that asymptotically q = i/t and hence G —> ^jqq, the boundary term in the improved action using (12.101) 
and (12.81) vanishes, 


y dciq 


- G 


= 0, 


o 


( 2 . 20 ) 
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thus making r[g c j] finite. It is also worth noting that ((2.1811 leads to 

dG \/2m x /Eq 2 + k 

dq q 2 

2 

and since asymptotically E \t f ^oo= 2 ^ 7 , one gets 


( 2 . 21 ) 


dG 

dq 


■P, 


q—t 00 


( 2 . 22 ) 


which is an explicit way of seeing dr[g c ;] = 0 . 

Next we analyze how the polymer quantization affects the above argument. For this we consider two polarizations 
of the polymer representation for the analogue model in the following sections. 


III. THE POLYMER MODEL WITH DISCRETE POSITION 


Polymer quantum mechanics is based on the idea of using the polymer representation of the Weyl algebra, a 
singular representation that does not obey the Stone-von Neumann theorem and hence is not equivalent to Schrodinger 
representation fl2l] . As is well known, the Weyl algebra is based on the Weyl relations 

KK = ^M 1 +A »2 1 (3-1) 

KK =Va 1 + a 2 , (3.2) 

%V X =e-^V x %, (3.3) 

with the reality conditions t/t = U_ fl and = V-\. In the polymer representation of this algebra, the corresponding 
Hilbert space, 'Hpoiy, possesses an uncountable orthonormal basis such that 

{lA v ) = <W, p, (3.4) 


where „ is a Kronecker delta. 

One can consider two polarizations of this polymer representation in which either the representation of [/ M or 
V\ on Hpoiy is not weakly continuous in its corresponding parameter, p or A. More precisely, by saying, e.g., the 
representation of V\ is not weakly continuous with respect to A, we mean lim.\_ > o(A t I^A|A < ) 7 ^ (mI^a=o|A 4 )- A similar 
criterion applies for a polarization in which the representation of U p is not weakly continuous with respect to p. Now, 
in a polarization where the representation of V x is not weakly continuous, the basic operators U^, V x act on the basis 
vectors |g) as 


Uy\q) =e* pq \ q), ( 3 . 5 ) 

V x \q)=\ q -X). ( 3 . 6 ) 

Since [/ M is weakly continuous in p in this polarization, one can write Up = e^ 9 . However, this is not the case for 
V\. Namely, since the representation of V x is not weakly continuous with respect to A, the momentum p cannot be 
represented on the Hilbert space T~L po iy in a well-defined manner as the generator of V x . Thus, although classically 
one may write V x = e^ Xp , this is not the case quantum mechanically and V x should be seen as an operator on its own 
and not as the exponentiation of the generator p. 

Furthermore, since we cannot take the limit A —► 0 due to the singularity of the representation, and once we consider 
A as a fixed free parameter of the theory, one can see from (13.61) that by starting from a certain q = qo, states get 
restricted to a (one-dimensional) lattice in q space where the wave functions 'k(g) = (q'l’F) have nonvanishing values 
only on the lattice points {q n \<ln = q 0 + n\, n € Z} for that qo; here we choose qo = 0. Thus we say q is discrete in this 
polarization and write the basis \q n ) as a countable one, labeled with n and the value of the momentum is restricted 
to — < p < We call this polarization with q discrete, q polarization. The corresponding Hilbert space T-L qo is 

only a superselected sector of TL po i y , such that "Hpoiy = ©o<g 0 <A ^90 ■ 

In another polarization, one in which is not weakly continuous, things are the other way around: while we 
can write V x = e^ Xp (and also classically C/ M = e^ M9 ), we may not write as an exponentiation of q, since q 
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cannot be represented on the Hilbert space 'Hpoiy in a well-defined manner as the generator of U p due to the criterion 
lim /i _ i> o(A|t/ A1 |A) ^ (A| [7^=01 A). In this case the basic operators U p: V\ act on the basis vectors |p) as 

U^p) =\p-p), (3.7) 

V x \p) = e ^ Xp \p). (3.8) 

Here, one can see from (ED that by starting from a certain p = po, the states are again restricted to a (one¬ 
dimensional) lattice in p space where the wave functions \H( p ) = (p(’I') have nonvanishing values only on the lattice 
points {p n \pn = po + n/Lt, bGZ} for that po ; here also we set po = 0. So we see that p is discrete in this polarization 
and write the basis \p n ) as a countable one, labeled with n whereas the value of position is restricted to — ^ < q < —. 
We call this polarization the p polarization, in which p is discrete. The Hilbert space T~L Po now is a superselected sector 
of 7/.poiy, such that "Hpoiy — CDo<p 0 </; Afpo ■ 

At this point a remark on notation is in order. In the rest of this work, whenever we use the q polarization, we 
adopt the notation \q n ) for the discrete position basis and | p) for the continuous “momentum” basis. On the other 
hand, for the case of the p polarization, the discrete momentum basis is written as \p n ) and the continuous position 
basis as | q). Additionally we emphasize that from now on we are going to work in the separable superselected Hilbert 
spaces 7~L qa or H po , and not in the full polymer Hilbert space. 

Let us first consider the q polarization in which 


q\q n ) 
V x \q n ) 

q\p) 

v x \p) 

Notice that (q n \p) = yj^e ~' qnp/h , /Ju dp\p)(p\ = 1, ai 
Now, classically the Euclidean action can be written e 


Qn\Qn) 5 

(3.9) 

\Qn - A), 

(3.10) 

-d P \p), 

l 

(3.11) 

ei Xp \p). 

(3.12) 


E„ex \qn)(q n \ = 1 (see Appendix®. 


Se = J dr (^p(J^j ~ H E ( q ,p)^j 

with 

He = - W(q). 

Using (13.101) . the kinetic term in this Hamiltonian can be represented as (see Appendix iBl) 


(3.13) 


(3.14) 


P 2 ^Pl = ^l( 2 -^-V-x)- (3-15) 

The potential -i? in this case can be represented using a regularization following Thiemann (3fi| 

-j==-V. x {xM, v x} 

= ^{v^’M + { v ^’MeT’ (3 - 16) 

where in the second line, we have chosen a specific synnnetrization. It is obvious that other types of orderings are also 
possible. The full Euclidean Hamiltonian (13.141) can be represented as (using the Dirac prescription {•, •} —> — i /n[-, •]; 
see Appendix [Dll 


He 


H 2 

2m\ 2 



- k 




VW\,Vx 


VW\,Vx 


V-, 


(3.17) 


One can act the above Hamiltonian on basis vectors \q n ) to get (Appendix iDl) 


He\ q n ) 


h 2 

2mA 2 


kh 4 

"F" 


Qn 


— A| — a/1'7™ + A|^ | Qn) ■ 


(^Wn) — \qn — A) — | q n + A)) 


(3.18) 
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A 4 W T 
kfi 4 



Figure 1. The form of potential after polymerization with A = 1 in q polarization. Note that for simplicity we have plotted 
q — 1| — sj\q + 1| j , a continuous graph, while it actually should be a discrete one. The maximum is finite and 


w 1 

JknF/JJ) 
here is equal to 4. 


Now the potential is not singular anymore at q = 0 as can be seen from (13.1811 above, and it is illustrated in Fig. [T] 
Also if the energy of the system is smaller than the peak of the potential and the particle’s initial position is such 
that 0 < go < q(W max ) where g(VF m ax) is the position corresponding to the potential peak, the semiclassical issue 
mentioned above is solved: since q remains finite at all times, then action if functionally differentiable (boundary 
terms vanish) and S[q c i] < oo. Therefore there is no need to add a counterterm to the action to be able to do a saddle 
point approximation. However, in a more general case, when the energy of the system is greater than the peak of the 
potential, we have 


^Loo ^ 00 ( 3J9 ) 

and due to (12.101) , we will still have the problem that the action will not necessarily be functionally differentiable and 
that even so, SE[qci] will not be finite. 

This analysis shows that this kind of polymerization, which does not bound the position but discretizes it, re¬ 
moves the singularity of the potential at the origin. However it does not solve the problems with the saddle point 
approximation, and we still need to add a Hamilton-Jacobi counterterm to get a well-defined action for this kind of 
approximat ion. 

One solution to both of the above problems (both functional differentiability of the action and/or its finiteness on 
classical solutions) comes from modifying the potential at infinity such that qt^oo < oo and consequently Sq = 0 at 
the (time) boundary. This is the subject of the next section. In the remaining part of this section, we investigate 
the effects of polymerization on the Hamilton-Jacobi counterterm in q polarization. To do this, we first derive the 
effective Euclidean Hamiltonian in q polarization using the path integral method. This turns out to be (see Appendix 

ED 


H 0 s 


2 h 2 
m\ 2 


sin 2 



-W T 



(3.20) 


where W T is the continuous counterpart of the potential derived using Thiemann’s regularization. A note about 
some subtleties is in order here: we could not directly argue that the form of the effective Hamiltonian will be (13.201) 
based on the form of the potential in (13.181) . This is so because if one computes the action of the full quantum 
Hamiltonian on \q n ) states, one gets the above effective potential but with discrete q n , and also the kinetic term will 

not be Ajj- sin 2 anymore since this is the form of the kinetic term when it acts on | p) states and not the |g n ) 

ones. On the other hand, if we use the | p) basis for the action of the full quantum Hamiltonian, we would not get 
the effective potential in (13.201) . In the first case (acting H on |g„)), one will get an additional f dp in the kinetic 
term (and also the potential will be discrete), while in the second case (acting H on |p}), there will be a in the 
potential term. The way to overcome these problems and get the effective action (13.201) is to use the path integral 
formulation and not something like H e g = (q n \H\q n ). What makes the path integral method useful is that one can 
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bring the extra out of the exponential of the action and turn it into an integral that appears in the measure of 
the path integral (see Appendix [Cl for more details). 

To find the explicit form of the polymer G po iy for this effective Hamiltonian, we need to solve (12.1211 . Considering 
the above Hamiltonian and noting that p = —, Eq. (12.121) becomes 



<9Gpoly 

dt 


(3.21) 


Again since the Hamiltonian does not depend explicitly on t, Eq. (12.121) implies that the derivative dG p°iy/dt must be 
a constant H = E, and an additive separation of variables in G po i y can be performed 


This turns (13.211) into 


or 


G P oi y (qS E, t ) = n(q, E) - Et. 


2 h 2 
toA 2 


Sill 



-E = 0 


(3.22) 

(3.23) 


on 

dq 




mkti 2 

2A 2 


- M - V\Q + A l) ^ • 


(3.24) 


As clearly seen from this, finding H in this case is much more involved due to the presence of sin -1 A way 

around this difficulty is to use a perturbative expansion of the right hand side of the above, around Y q = 0. This 
yields 


an = Vy mVE? + * V / 2 „„ 3 ktfV 2 Vm\ r V 

dq q 12 H 2 q 3 ^ V ’ yjEq 2 + k J \h 4 


(3.25) 


The leading continuum order in (I3.25D . in which A = 0, matches exactly the corresponding classical counterpart in 
(12.171) . Integrating (13.2511 and then substituting the result into (13.221) yields 


G P oi y = — Et + \f2yfni\J Eq 2 + k + \f2\fmk tan 


-l 


Eq 2 


+ 


h 2 


V%\fmyj Eq 2 + k 

24Y 


(to (2 Eq 2 — fc) — 3S 2 ) + 


\f2sfm 

8 y/k 


E(+mk + h 2 ) tan 


Eq 2 + k 


(3.26) 



(3.27) 


Again we see that the purely classical term with A = 0 matches exactly to the classical Hamilton’s principal function 
(12.181) while there are also several types of corrections due to polymer quantization. 

Since the classical part of the polymerized G po i y in (13.271) is exactly the same as the classical nonpolymerized 
case, and since we have seen that with purely classical G, the action is functionally differentiable and its value on 
classical solutions is finite, we conclude that the polymerized improved action Tpoiyfgci] with the above counterterm 
has also the same nice properties and thus makes it possible to proceed with the saddle point approximation if desired. 
But in addition to that, the counterterm (13.271) has additional terms proportional to the powers of “quantum lattice 
parameter” A. So, it is reasonable to expect that this, together with the change of the form of the bulk action due 
to polymer quantization, will change the thermodynamical properties of the system in case of, e.g., a black hole, 
and therefore it is interesting to see what are the implications of such semiclassical polymer modifications. Next we 
consider the other polymer polarization. 


IV. THE POLYMER MODEL WITH DISCRETE MOMENTUM 

In this section we consider the p polarization in which 


U»\Pn) =\Pn - lA, 

(4.1) 

P\Pn) =Pn\Pn), 

(4.2) 

UM = e ^ 9 k)> 

(4.3) 

p\q) =-^d q \q). 

(4.4) 
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- 2 

The representation of the kinetic term in (13.141) in this polarization is very simple; in fact it is just 4^-. The problem 
here is how to represent the potential. It is not clear how Thiemann’s regularization can be used in this case. The 
reason is that generally this regularization is used to represent a variable that is discrete and not bounded. In the 
present case though we have to define a replacement for the inverse of q\ a problematic task since even q is not well 
defined on "H Po but only U lt . Even more, we may consider finding functions _F(t/ p ) and G{p) such that classically 

-L = {F(f7 p ) ,G(p)} m , n,m> 0, (4.5) 

q n 

so that F(Un) and G{p) admit a simple representation on Hilbert space. The first part, i.e., finding classical functions 
F(l Jfj) and G(p) fulfilling (14.51) may not be very hard, and several options may be available such as 

A {v'l - |,p} = F {vWM.p} = <«> 


However, the second part, the ability to represent the functions F(U p ) and G(p) on Hilbert space, is the hard part 
as can be seen from the above example. Presently we have not found satisfactory functions that can be represented 
on Hpoiy for which Thiemann’s regularization can be done. In spite of this difficulty it is still possible to introduce a 
formal inverse squared position operator in this polarization and hence its semiclassical approximation. 

We saw in the previous subsection that a solution to both issues of functional differentiability and finiteness of the 
action is likely to come from bounding q to finite values, and this may result from modifying the potential term in 
He- Considering this, we can argue that since in this polarization, q is bounded, a representation W(U^) is expected 
to lead to a bounded potential and thus would eliminate the need to add a Hamilton-Jacobi boundary counterterm to 
the action to enable one to make a well-defined saddle point approximation. Now we provide a scheme that shows how 
the polymer quantization will change the action in such a way that no counterterm is needed to have a well-defined 
saddle point approximation. 

Our proposed scheme for the effective potential is based on the observation that q 2 in this polarization gets replaced 
by an operator that is represented as (see Appendix dBl) ') 

Q 2 ^ (2 - % - U-,) ■ (4.7) 

Its action on a | q) basis is 

This is computed using (14.31) (see Appendix (0). Based on this, our proposed replacement of the inverse squared 
operator is 



Using path integral formulation, the effective action turns out to be (Appendix [Cl) 


S e ff 





(4.10) 


Eq. (14.101) suggests that the effective form of the classical potential W(q) = 4)- in this scheme has an effective form 


W h 



(4.11) 


Note that here also there are subtleties similar to those explained in the previous section related to the important 
role of the path integral method to derive (14.101) . For example, acting quantum Hamiltonian on | q) yields the above 
effective potential while the kinetic term will not be p 2 /2m anymore since this is the form of the kinetic term when 
it acts on | p) not on | q) . Furthermore, since in this polarization our momentum states are actually discrete, i.e., 

| p n ), p n = np,n £ Z, even if we act the kinetic term on them instead of | q), we will get a discrete result and not a 
continuous one. Thus again the path integral method saves the day and makes it possible to get rid of the additional 
V., and the discreteness of the momentum in the effective action. Details are explained in Appendix [C] 
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w h 



Figure 2. The form of potential W h for our model with //,= !. 


The plot of the potential (Id.Ill) is shown in Fig. [2j Obviously the potential is still singular at the origin, q = 0. 
But the good news is that the position is bounded at the time boundary, i.e., q c i(t —> oo) < oo. This means that 
since now q is finite at the boundary, Sq = 0 on the boundary, and thus the action is welldefined for a saddle point 
approximation without the need to add any other terms, such as the Hamilton-Jacobi counterterm. This polarization, 
like the previous one, will also most probably lead to modifications to the thermodynamics of the system due to 
the modification of the bulk action due to polymer quantization as well as the absence of any counterterm. This is 
particularly interesting if the system under study is a black hole. 


V. DISCUSSION 

In this work we have studied the issue of ill-defined saddle point semiclassical approximation that occurs for some 
systems including dilatonic black holes and some mathematically analog mechanical models. This problem arises due 
to the fact that the action is not functionally differentiable for all the variations of the fields, compatible with the 
boundary and falloff conditions and even if so, the value of the action on classical solutions does not remain finite for 
all of these conditions. This issue is rather important since one of the main methods of deriving thermodynamical 
properties of such models is through the saddle point approximation to the path integral that can also be interpreted 
as the partition function of the system. The common solution to this problem is to add a boundary counterterm to 
the action which is a solution to the corresponding Hamilton-Jacobi equation. A very interesting observation is that 
only with this term can one get the correct thermodynamics for certain black holes. 

We undertake to seek an alternative method to attack this issue which is through polymer quantization of the 
model. The effects of polymerization may lead to two outcomes: either it will spare us altogether from adding any 
additional term to the action and the action is already well defined for saddle point approximation after polymer 
modifications, or it may not remove the necessity of adding a counterterm, but it will modify it. In both cases, it 
is very likely that the process of polymer quantization will change the thermodynamics, either due to the fact the 
counterterm is not necessary or because it will be modified. 

In this work we restrict our study to a simpler analog model featuring the aforementioned issues. This helps us see 
the effects, problems, and results of this method more clearly and provides interesting insights about the procedure. 
The analog model we study is a single particle in an inverse square potential. 

We first show how this system is modified under different polarizations of polymer quantization. It turns out that in 
the polarization where q is discrete, we still will need to add the Hamilton-Jacobi counterterm to the action, but the 
advantage is that the potential can be rather easily represented on a Hilbert space using Thiemann’s regularization 
f36i |. This is the case in which the Hamilton-Jacobi counterterm is modified by polymerization. We then proceed to 
compute the effective action and thus derive the effective Hamiltonian using path integral formulation. Then we derive 
the associated Hamilton-Jacobi counterterm and show that the classical terms of the polymerized case match exactly 
the nonpolymerized case while there are corrections to this counterterm that come from the polymer quantization. 
This supports our claim that polymer quantization will change the thermodynamics of the system due to the polymer 
modifications of the bulk action and the Hamilton-Jacobi counterterm. 

On the other hand, in the polarization where p is discrete, we argue that in general there should be no need to 
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add a counterterm to the action since the variable q is bounded due to polymer effects; thus the action remains 
functionally differentiable and finite on classical solutions in accordance with all the allowed variations. However, 
since the representation of the potential in this case is not so straightforward and could not be done here, in order to 
get concrete results and calculations, we proposed an effective form that replaces the classical potential W = -^ based 

on the analysis of the polymer operator q^. Using this, we show that the effective action is indeed finite evaluated on 
the effective solution and is functionally differentiable without the need to add a counterterm. This means that due to 
the effects of polymer quantization in p polarization, the system is already well defined for saddle point approximation. 

Further developments along the lines of the present work include the following. It is pretty evident that a similar 
analysis for the case of polymer black hole systems may lead to the change of thermodynamics due to polymer 
modifications of the bulk action. It will be interesting to check whether any counterterm is required in such a case. 
Also, recently, much work has been done to unveil singularity avoidance in loop quantum cosmology 0 as well as in 
black holes [37H39| by alluding to effective models like the ones we have studied in the present work; hence it would 
be interesting to explore possible physical consequences for a semiclassical approximation of these systems using the 
path integral framework. 
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Appendix A: Proof of (13.121) and (14.31) 


Let us consider the polarization (g, V\) where eigenvalues of the operator q are discrete. We would like to derive 
(13.121) from (13.91) and (13.101) . Using the completeness relation 


and the form of 


one can write 


1 = 


\ /■"'‘A 

— / dp\p){p\ 

£TTtl J — ttH/x 




Vx\ p) = V x \q n )(q n \p) 

n =—oo 

\ /•"A 00 / \ 

= 2£ llWkn-A)^/—e* 

r / x oo . 

/ dp' J2 \p')e~* pqn e* pX e* pqn 
J —’’’/A _ 


PQn 


X_ 

2-7T 


=A 


/ 


-7t/A 

/ n =—oo 

t A 


dp' lim 

-*A 


— y 

2tt ^ 


(p-p‘) 


AP' X In' 


IA 


<5 n(p-p') 


(Al) 


(A2) 


(A3) 


where 6 n (p — p') is a Dirac delta sequence and the limit will make it into a formal Dirac delta series representation 
S(p — p'). Using the above and the fact that for a function f{p) with a period 2v /\ one can write 


[ dp'8{p - p')f{p') = \f(p), 

J -*a A 


(A4) 
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we finally arrive at 

V\\p) =e% pX \p). (A5) 

Using the same lines of argument but now in the polarization (U^,p) where eigenvalues of the operator p are discrete, 
one can derive (14.311 from (14.11) and (14.211 . 


Appendix B: Action of q^\q) and p\\p) 


Let us choose one of the polarizations, say the one in which p is discrete. Since as mentioned before, the q operator 
does not exist in this polarization, we first need to construct an analog of the q 2 operator that we call q 2 and then 
find its action on the desired basis. To do this, we note that since classically U^ = e M t n , then one can write 


2 - 


pV 

K 2 ’ 


h 


(Bl) 


Using this, one can isolate q 2 and represent its singular counterpart q 2 (singular here means the limit p —> 0 does not 
exist) on the Hilbert space as 

H 2 


qp p 2 


2 -Up- U-„ 


)■ 


and using the results of Appendix [A] or equivalently using (14.31) , one can write 

h 2 r. \. h 2 

,2 




2-U,-U.,)\q)=^ 

h 2 


-)»-£(■ 


2 — e hM — e ~hL 


W) 


hr / 

=—z (2 — 2 cos 
V 


m '<> 


So we conclude that 


ql\q) = (sin 2 (§)) |g). 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 


On the other hand, if we choose to work in the polarizations in which q is discrete, the operator p\ can be defined 
and represented using the same method, i.e., 

h 2 


vl = 


— ( 

A 2 V 


2 - Ua - U_ 


0- 


(B7) 


Its action on the \p) basis can then be computed in the same manner as above for q 2 , and by applying the results of 
Appendix [A] or equivalently Eq. (13.121) one gets 

Pl\p) = % (sin 2 (^)) Ip)- (B8) 


Appendix C: Effective polymerized actions 


In this section we show how to obtain the effective action of the models described in Secs. IIIIl and lTVl The calculations 
are done in Newtonian form but changing to the Euclidean form is straightforward. In steps that are different for 
each polarization, we will make comments and make clear what equation corresponds to which polarization. The 
transition amplitude in each case is given by 


((?/; tf I ti) 

(Qf: tf | qi 1 ti ) 


lim 

N—too 


lim 

N—> 00 


( N—l 


N 


nz n | Qn— l?^n— l) I ^ 


\ n— 1 n 
/N—l n l 


\n—l 


N 


n 


dq n 


| Qn—lit n —l) J , 


(Discrete q), 

(Discrete p). 


(Cl) 

(C2) 
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where in the first equation, we have a sum over polymer lattice q points due to discreteness of q and in the second 
equation, the limits of the integral reflect the bounds on the continuous variable q. Note that in the first equation, 
the subscript is due to the genuine polymer lattice discreteness of < 7 , while in the second one, it is due to partitioning 
of the path integral. The amplitude above is then divided into partitions for which we have 


{q n ,t n | q n -i,t n -i} ={q n | e ^ poly | q n - 1 ) 


l€ . 


| 1 -f^poly | Qn— l) 

/ , x (2-V x -V-x 

-(Qn I qn- 1 } h (qn | ^ ^ 


qn— 1 } ^ (qn | W | qn— l)j 


(C3) 


where H po \ y is the polymer quantum Hamiltonian operator and we have used (1B7I) for the representation of the 
kinetic term. W is the potential operator that can be represented by Thiemann’s regularization as in (13.1711 in case 
q is discrete or can be the potential proposed in Sec. m when p is discrete. We then insert one of the following 


identities: 


1 = ^2\Pn){Pn |, (Discrete p ), 

n 

-rrh 

A f— 

2ttH J -nh dpn \ pn }& n \' (Discrete q), 


1 = 


(C4) 

(C5) 


in front of the kinetic and potential terms in (IC3I) . corresponding to the polarization we are working in, to obtain 
either 


7r n 

A f-r 

{Qm^n | Qn—li^n—l) ~ 2irh / ^ Pn{,Qn \ Pn){Pn \ Qn— l) 


i-sr 


2 h 2 


sin 2 I ^'1 +w .: j 


h \m\ 2 \2h 


(Discrete q ) (C 6 ) 


or 


{qn^n | Qn— lj^n—l) | Pn){Pn | q n - 1) 


1 - — 

fPL + W h\ 

h 

{2m n ) 


(Discrete p). 


(C7) 


Here 


Wn (a/I qn - A| - \/\q n + A|^ , 

W h = ^ csc2 (^\ 
n 4/i 2 V 2 h) 


(C 8 ) 

(C9) 


are Thiemann-regularized and heuristic polymer potentials, respectively. Note that again there are two types of 
subscripts associated with two types of discreteness. One is related to partitioning the full transition amplitude [as 
in and the other is related to the genuine polymer lattice discretization [as in (|C 8 [) j. Using 


<» i-”•>='/ 




(Discrete q), 
(Discrete p), 


and also the full expansion of the exponential, one gets for the above partition transition amplitudes 

2 


{Qm^n | Qn—l't^n—l) — ^ 2ttH ^ 

('Qni tn | Qn— 1? tn— l) = ^ ^ e 




T 


[pn <a "-; n - l) - (^T 'In 2 (^D+W^)] 

(Discrete p). 


dp n e^y p 

— -(lt+<)] 


(Discrete q), 


(CIO) 

(Cll) 


(C12) 

(C13) 
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Using these in corresponding full transition amplitudes (IC1I) and (IC2I) yields 


/N -1 


( qf ,tf\ qi ,u )= N ^^ 0 IIE 


(A) 


2 N 

n 


dpn e 


Etif[p 


(.Qn-Qn-l) ^ 2h 2 ( Xp 


n= 1 17 __ r 


f AT—1 


{Qf,tf\qi,ti) = lim 

N — >-oo,e—>-0 


n 

v, n—1 


dq„ 


P 

27tH 


N 


iie 

n=l n / 




- s+w: 


)+ w n)] (Di scr ete q), 

(C14) 

(Discrete p). (C15) 


In both cases above, we have aj] in front of the exponential while we should have an integral instead. In other 
words, if we wish to have a normal path integral, we need to pass from a discrete variable, say g„, to the continuous 
one q (and the same for p n ). To achieve this we use the identity 0 


p2ll pOO pOO 

E / dqf{q,p n )e lPnq = / dq dpf(p,q) 
„ J 0 J — oo J— co 


(C16) 


valid for ,f(q,p) which are periodic in q, and a similar expression also holds for the case with discrete p n . Let us 
consider first the expression (1C 151) . Using the above identity, one gets for (IC15I) 


(q f ,t f \qi,ti) = lim 

N — »oo,e— 


f N—l 

n 

\n=l * 


dq n 


P 

27 vh 


N 

n 

n—l 1 


dp e 


ie v'IV |„ Qn — 1 / y n , >ir 

~h En = l IP™ -c-( 


*)] 


(Discrete p), (C17) 


which is the polymer Feynman formula. Taking the limits in (1C 171) above yields 


(qf,t f \qi,ti) = ^ (il / dq ”) (il/ dp\e hSdt \- pq ( 2m+W )] (Discrete p), (C18) 


where the effective action can be read off to be 


S'eff = J dt [pq - H c ff ] = J dt 


pq 


P prk 2 

2 ^ + W° SC 


v 2 h) 


(Discrete p). 


The case with discrete q is a bit more delicate because of the inherent discreteness that precludes the limit 

Qn Qn—1 . 

->pq- 


lim p n - 

e—>0 1 


(C19) 


(C20) 


To deal with this, one uses a method similar to Leibniz rule for discrete variables. One can rewrite the sum that 
appears in the kinetic term as 


N 

E Pn 
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qn qn—i l 
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N 


N 


=- (E pn?n ~ E^ 5 "- 1 

n—l / 

N—l 
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(C21) 


which is similar to expressing pq = —qp+ (pq) where - ( PNqN — Piqo) plays the rule of the “boundary term” (pq)' in 
a discrete sense. Using this, expression (1C14I) becomes 
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Now, we can take the limit such that the term lim c _>o ^ p " +1 pDgr.+pwgiv pi<?o becomes — qp + (pq), which in the 
continuous limit can be rewritten as pq. Before taking this limit, we use an identity similar to (IC16I) but for discrete 
q , 

p2lV pOO poo 

E / dpf(q n ,p)e i ^= dq dpf(p, q)e iqp , (C23) 

„ J 0 J— oo J— oo 


to turn into an integral. Then by taking the limits, one gets 


{qf,t f \qi,U) 





sin 2 (M) +w ^)] 


(Discrete q) (C24) 


Note that the potential does not have a discrete subscript anymore since we turned q n into a continuous variable q 
by using (IC23I) . Thus the effective action can be read off to be 


S'eff = J dt [pq - H eS \ = J 


dt 


PQ' 


2 h 2 
mX 2 


sin 




(Discrete q). 

(C25) 


Appendix D: Action of He in (q,V\) polarization 
1. In | q„) basis 


The action of the Hamiltonian (13.171) on the basis \q n ) in the polarization (I3.9D - (I3.12I) can be computed as follows. 
The kinetic term acts like 


(2-Vx-V-x\. , h 2 . . . . . 

2m ~~ 2m ( A* ) |9n) ~ 2^? (2kn) Wn A) l<Zn + A)) ’ 


(Dl) 


where the operator p 2 x is defined in (IB7I) . The potential W in (13.171) has two terms for each of which we have 
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Thus the full expression for the action of the potential in this basis becomes 
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(D2) 


(D3) 


(D4) 


2. In |p) basis 

One can act the Hamiltonian (13.171) on basis |p) where 

V\\p) = e% Xp \p). 


Then the kinetic term in (13.171) turns out to be 
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where we have used dm ) and (IB8I) . To find the action of the potential term, we first note that 
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Then using the results of the previous subsection, we get for the action of the potential in this case 
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